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Abstract 

In a previous article, an 'invariant method' to calculate monomial 
integrals over the U (n) group was introduced. In this paper, we study the 
more traditional group-theoretical method, and compare its strengths and 
weaknesses with those of the invariant method. As a result, we are able to 
introduce a 'hybrid method' which combines the respective strengths of 
the other two methods. There are many examples in the paper illustrating 
how each of these methods works. 



1 Introduction 

This article deals with the calculation of integrals of the form 

(du)u* in ---u; pJ u klh ---u kqlq (i) 

over the U(n) group, where (dll) is the invariant Haar measure normalized 
to J(dU) = 1, and Uij is a U(n) matrix element, with being its complex 
conjugate. 

These integrals and their generating functions are useful in many areas of 
physics, including two-dimensional quantum gravity £Q, QCD, matrix models, 
and statistical and condensed-matter problems of various sorts They are 
also needed in the parton saturation problem at small Feynman-cc [3]. 

The integral JIJ depends on the indices I = {i\ ■ ■ ■ i p }, J = {j\ ■ ■ -j p }, K = 
{fci • • • k q }, and L = - ■ ■ l q }, so it will be denoted as (IJ\KL): 

(IJ\KL) = f(dU) UljU KL , (2) 



where Ujj = Ylo=i U* a ,- j an d similarly for Ukl- Since the matrix elements 
commute, Ufj — Uj pJp where Ip = {ip(x) • ■ 'ip(p)} is obtained from / by a 
permutation P G S p of its p indices. Hence 

(IJ\KL) = {I P Jp\K T L T ) (3) 
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for any P G S p and T 6 S q . 

The integral is nonzero only when p — q, a number which will be referred 
to as the degree of the integral. Without loss of generality, it turns out that we 
may assume K = I and L to be a permutation of J, namely, L — Jq for some 
Q G S p . The value of the integral depends on what the index sets I, J, and 
what the element Q are, so even for a given p, there are many distinct cases. 
The best way to distinguish them is to represent each integral by a diagram in 
a way to be explained in the next section. 

Integral JTJ has been computed using a graphical technique 0]. It can also 
be obtained using the Itzykson-Zuber formula |S] as a generating function, or 
directly from group theory [S] using the Frobenius formula |7j. We shall refer 
to this last method as the group-theoretical method, or GTM for short. In the 
GTM, a general formula is available to compute (|TJ. It involves a triple sum 
over an expression containing characters of the symmetric group S p , as well as 
the dimensions of irreducible representations of S p and of the unitary group 
U(n). One of the sums is taken over all the relevant irreducible representations, 
and the others are taken over the symmetry groups of the index sets / and J. 
These sums could be long and tedious for a large p, and for most symmetry 
groups. 

A different way to calculate was introduced in a recent paper [8]. This 
method relies only on the unitary nature of the matrix elements, and the in- 
variance of the Haar measure. In particular, no knowledge of group theory is 
necessary. The invariance of the Haar measure, as well as the off-diagonal uni- 
tarity relation, are used to derive relations between integrals of the same degree. 
The diagonal unitarity relation connects integrals of degree p with ones of degree 
p — 1. Through a chain of these relations, the desired integral is finally related 
to the basic integral of degree 0, which is J dU = 1. The desired integral is then 
solved from this chain of relations. We have called this method the invariant 
method, or IM for short. 

The purpose of this paper is to compare the pros and cons of the GTM 
with the IM. In order to do so we must first study and understand better the 
nature of the GTM. Armed with this comparison, we will be able to design a 
new hybrid method which combines the strengths of these two other methods. 

The IM is reviewed in Sec. 2. It is used to derive a new 'double-fan' relation 
needed in the example for the hybrid method. 

The GTM is reviewed and studied in Sec. 3. It can be used to derive simple 
relations, but the relations derived in this way are nowhere as powerful as those 
derived with the IM. The group-theoretical formula can be used to calculate 
any integral, but generally that is tedious and has to be done integral by in- 
tegral. However, for integrals whose two symmetry groups are disjoint, or one 
is contained in the other, systematics emerge to make the calculation simpler. 
Several of these 'orderly' integrals are studied in Sec. 3. 

A comparison of the strengths and weaknesses of the two methods is to be 
found in Sec. 4. Armed with an understanding of their relative merits, we design 
a 'hybrid method' in Sec. 5 to take advantage of their respective strengths. This 
method is illustrated by the calculation of a class of 'double-fan' integrals. There 



2 



are also four Appendices showing details of various calculations. 

2 The Invariant Method (IM) 
2.1 A Brief Review 

The invariant method presented in a previous paper [5] can be used to calculate 
the integral QJ. The method exploits the unitarity of the U(n) group elements, 

n n 

and the invariance of the Haar measure, in the form 

J(dU)f(U,U*) = f(dU)f(U*,U) = J(dU)f(U T ,U* T ) 

= J{dU)f(VU,V*U*)=J{dU)f(UV,U*V*), (5) 

for any function / and any V E U(n). 

The unitarity relation Q relates integrals of the same degree if i ^ I, and it 
relates integrals of degree p to integrals of degree p—1 if i = I. Other relations 
between integrals of the same degree can be obtained from (JSJ, by suitable 
choices of V. Here are some of them discussed in the previous paper 

1. Using Vij = e^Sij, it follows that p must be equal to q in order to avoid 
the vanishing of integral J5J). The number p will be called the degree of 
the integral. 

2. Using Vij = e^Sij, it follows that J2J) is nonzero only when K = Im and 
L = Jr for some M, R € S p . Using ©, and denoting RM^ 1 by Q, only 
integrals of the type 

(IJ\IJ Q ) = j \dU)U* u U IjQ (6) 

are nonzero, so from now on we need to consider only integrals of this 
type. 

An integral with Jq — J will be called a direct integral. For such integrals 
we may always choose Q — e, the identity permutation. Otherwise, the 
integral is an exchange integral. 

Integrals are represented diagrammatically as follows. Each distinct value 
in the index set I is represented by a dot on the left (L-dot), and each 
distinct value of the index set J or Jq is represented by a dot on the right 
(R-dot). The factor U*j is shown as a thin solid line between the L-dot 
i and the R-dot j, and the factor E/y is shown as a dotted line between 
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(a) (b) 

Figure 1: Examples of U(n) integral diagrams, (a) The unique exchange integral 
for p — 2; (b) a Z-integral with arbitrary multiplicities mi, m-z, and 7713. 



these two dots. The factor U^Uij is represented by a thick line, or more 
generally, the factor U*^jUh is represented by a thick line with a pair of 
numbers (m, n) written beside it. If m — n, then only a single number m 
is written. The numbers (m, n) or m will be known as the multiplicities 
of the line. See Fig. ^ for an illustration. 

3. With V chosen to be a permutation matrix of n objects, it follows that 

(IJ\KL) = (I'J\K'L) = (IJ'\KL'), 

where I' is obtained from I by a reassignment of the values of its indices, 
e.g. if / = (334), then I' may be (558). K' is obtained from K by the 
same reassignment, and similarly for J' and L' . As a result, there is no 
need to know the values of the indices of the L-dots, nor the R-dots. This 
is why the dots in Fig. are not labelled. 

4. As a consequence of the first two equalities in J5|, an integral remains the 
same under the interchange of the solid lines with the dotted lines, or the 
L-dots with the R-dots. 

5. Using V = R(ab), the rotation matrix in the (a, b) plane, a 'spin-off re- 
lation' is obtained. Consider a R-dot imbedded in an arbitrary integral 
Mo = (IJ\IJq), with d pairs of solid-dotted lines attached to the dot. Now 
spin off e pairs of these lines to create a new R-dot and a new integral. 
There are many ways to choose the e pair of lines, each possibly corre- 
sponds to a different integral. Let M e be the sum of all these integrals. 
Then the quantities Mo and M e are related by the spin-off relation 

M e = M (f), (7) 
where f | is the binomial coefficient. 

The relation is local in that it is independent of the the structure of the 
rest of the diagram. The same relation can also be used to spin off a L-dot. 
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(a) (b) 

Figure 2: The fan diagrams shown here can be a part of a larger diagram. In 
that case, there may be many more dots and many more lines in the complete 
diagram, provided none of the additional lines land on the R-dots shown, (a) A 
closed fan; (b) a partially opened fan. If all m» = 1, then it is said to be a fully 
opened fan, or just an opened fan. 



In what follows, we summarize two general results obtained in the previous 
paper [Sj using the IM. 

2.1.1 The fan relation 

The fan relation 

(dU)A\U acl \ 2m ^--\U aCt \^ = J " UdU)A\U ac \ 2m (8) 

J \2^,j=i m j)- J 

relates integrals of the same m = Y^j=i m h where A is an arbitrary prod- 
uct of matrix elements of U and U* whose column indices are different from 
c\ , C2 , . . . , Ct ■ The column index c on the right could be taken to be one of the 

Cj's. 

Diagrammatically, the integral on the right of (JSJ is shown in Fig.[2a), and 
the integral on the left is shown in Fig. 0{b). The additional lines and dots 
corresponding to the factor A are not shown, because they do not affect the 
spin-off relation. We shall refer to Fig. Efa) as a closed fan, and Fig. Efb) as a 
partially opened fan. If every mj = 1, then it will be said to be a fully opened 
fan, or simply an opened fan. 

In particular, a closed fan integral is m! times an opened fan integral. In 
fact, this relation between the two types of fans immediately gives rise to the 
relation (JHJ) between a closed fan and a partially opened fan. To see it, note 
that each branch of a partially opened fan is itself a closed fan. By opening up 
all of them, we get the fully opened fan integral, multiplied by a multiplicity 
factor Y[j m j! from all the branches. Thus a closed fan is ml/ FJ . rrij\ times a 
partially opened fan, as given by ||5J. 
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(a) 



(b) 



Figure 3: Double- fan diagrams, (a) A closed diagram; (b) a fully opened dia- 
gram. 



2.1.2 The Z-integral and the fan integral 

In we also obtained a general formula for integrals of the type shown in 
Fig.^b), for arbitrary nonnegative integers mi,m2, and m^. We call that the 
'Z-formula' because of the shape of the graph. It is 



(9) 



Z(m, 1: m 2 ,m 3 ) = J (dU)\U lJ \ 2m '\U tl \ 2m2 \U kl \ 2m3 

mil m 2 \ 7713! (n — 2)l(n — 1)! (n + mi + m% — 2)! 
(n + mi — 2)! (n + 7713 — 2)! (?i + mi + m 2 + m^ — 1)! 

In the special case where m 2 — m 3 — 0, this becomes 

F(m) ^ Z(m, 0, 0) = .^"^l' (10) 
(n + m — lj! 

which is the integral for Fig. EJa) when there are no additional dots or lines 
around. 



2.2 Double-Fan Relation 

The fan relation can be generalized to a double-fan relation, connecting the 
closed 'double-fan' diagram of Fig.^a) with a (fully) opened double-fan diagram 
such as Fig. EHb). As in Fig. [21 there may be additional dots and lines in the 
integral, but none of them may end up on the R-dots shown. From that relation, 
we can also deduce relations between a closed double-fan and a partially opened 
double-fan, as done in the single-fan case. 

The double-fan relation is considerably more complicated than the single-fan 
relation JSJ), because there are many more double-fan graphs. Each R-dot of a 
(fully) opened (double-fan) graph such as Fig. EJb) falls into one of four basic 
patterns: [A a ], [At,], [B a ], and shown in Fig.0] If the solid and dotted lines 
end up on the same L-dot, the pattern is a [£?]; otherwise it is an [A]. The 
subscripts a and b tell us which L-dot the solid line emerges from. 



G 



[AJ [Ab\ lB a ] [B b ] 

Figure 4: The four basic patterns for the R-dots of a fully opened double-fan 
graph. 



Suppose there are at number of [Ai] and f3i number of [Bi] patterns in a 
(fully) opened (double-fan) graph. Then there are m a solid and n a dotted lines 
emerging from the L-dot a, and m b solid and n b dotted lines emerging from the 
L-dot b, where 

m a = a a +f3 a , n a = a b + (3 a , 

m b = a b + f3 b , n b = a a +(3 b . (11) 

The total number of R-dots in the opened graph is N = m a + m b — n a + n b . 

When these TV R-dots are merged together, we get the closed (double-fan) 
graph depicted in Fig.EJa), which will be denoted by [(m a n a )(m b n b )]. If N = 1, 
this is just one of the four basic patterns discussed before. If TV > 1, we will 
call it a compound pattern. 

From (|TTj) . we see that if we replace ai and /J; by a- = +£ and = — £, 
with any integral £ which keeps a[ and (3[ nonnegative, then we get the same 
closed graph by collapsing this new opened graph. Conversely, it will be shown 
in Appendix IXI that the closed graph [(m a n a )(m b n b )} can be spinned off into 
a sum of several opened graphs, one for each (a^Pl). The double- fan relation 
expressing that quantitatively is 

{(m a n a )(m b n b )] = ^(a'^'M) [A a ]<[A b ]<[B a f«[B b f» 

= Y,KAa + a' b A b + [3> a B a + P b B b }, (12) 

where 

vMM - (13) 

and the sum is over all solutions (a' a , a' b ,(3' a ,(3' b ) of 111 II) . 

The double fan becomes a single fan if the L-dot b is not connected, namely, 
if m b = n b = and m a = n a = m. In that case l|12|) becomes 

[mB a ]=m\ [B a ] m , 

which is just (JHJ (when all mi = 1) in another notation. 
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3 Group-Theoretical Method (GTM) 



3.1 A Brief Review 

Using group theory, the integral © can be turned into a multiple sum jSJ |S] . 
In the notation used in Appendix A of jo], the formula is 

(u\Uq)= E E E7ifx X/(5Qi?) ' (14) 

Reg, see./ / W-) a f 

where p is the degree of the integral. The symbols Qi and Qj represent the 
symmetry groups of the row and column index sets. More precisely, Qx = 
{P G S P \P(X) = X}, with X being either I or J. The irreducible represen- 
tations of the symmetric and unitary groups are both labelled by a sequence 
/ = {fx, /2, ■ • ■ , f P ), with fx > h > ■ ■ ■ > fp > 0- Xf(P) is the character of 
P G S p in the irreducible representation with signature /. The dimension of the 
irreducible representation / is given by df = Xf( e ) f° r S p , and by df for U(n). 
A formula for df is: 

*f> ;.)- w ' + g-y' + '"- 3 > r - / - > . < 15 > 

— 1, n — 2, . . . , 0) 

where D(x\, . . . , x n ) is the Vandermonde determinant given by Il^fcC 2 -* ~ ^fe)- 
Since Xfi.9) depends only on the class c g that g belongs to, we may write it 
as Xf( c g)- With this notation, eq. (fT4*jl can be re-written as 

(/j|/j Q > = E^M£M, (i6) 

c 

where the sum is taken over all classes c of S Pi 

n w = E E S ^Q R G c ) 

HeSi seg., 

is the number of elements of the type SQR in the class c, and 

^]=E7rfx X ' (c) - (18) 

It is not difficult to see that Q is not unique, because Q' — SQ = QT for 
any S G Qj and any T £ 5j Q is another possible Q. It does not matter which 
Q we pick in (|17|) . That equation can also be written as 

= E E 5 ^ TR e c )- ( 19 ) 

Reg, Teg., Q 

It is straightforward but generally very tedious to compute N[c], because 
we need to calculate the product QTR for every T G Gj q , every R & Gi, and 
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determine what class c the product belongs to. Then we have to count up all 
the products that are in a given class c to get N[c\. However, the task becomes 
considerably more manageable if either Qj and Qj Q are disjoint, or if one is 
contained in the other. We shall refer to integrals with those properties as 
orderly. Further simplification occurs for direct integrals, because in that case 
Q can always be chosen to be the identity e, so the triple product is reduced to 
a double product TR. 

The calculation of £[c] in l|18|) is simpler than the calculation of N[c], but 
still we know of no closed form of it valid for every class c and every symmetric 
group S p . The best we can do is to compute them case by case. Results are 
given in Sec. 3.3. Each £[c] is actually an orderly integral with Qi = Qj = e, to 
be referred to as a primitive integral. 

Other integrals can be computed in terms of the primitive integrals, if N[c] 
is known. We shall discuss two orderly integrals for which N[c] can easily be 
obtained. In Sec. 3.4, we discuss the stack integrals, which are direct integrals 
with Qj = Qj. In Sec. 3.5, we discuss the fully opened double-fan integrals of 



Relations between orderly integrals may be obtained without knowing the 
explicit values of £[c], if their 7V[c]'s are related in a simple way. This is the case 
for the single-fan relation, and the double-fan relation with n a = = and 
m a = n b = m. They will be discussed in Sec. 3.2. However, general double-fan 
integrals are not orderly, so we cannot obtain the general double-fan relation 
by the group-theoretical method, at least not in the present way. We will also 
show that the closed (single-)fan integral can also be computed without explicitly 
knowing what £[c] are. This is one of the very few cases where integrals can be 
obtained group-theoretically without explicitly knowing £[c]. 

That leaves the non-orderly integrals, for which each term of the summand 
in 119|l has to be calculated separately to get N[c]. The first non-orderly integral 
occurs in degree p = 3. In Sec. 3.6, we shall show how to calculate some of the 
p = 3 and p = 4 non-orderly integrals. 

3.2 Single- Fan and Simple Double-Fan Relations 

The single-fan integrals are orderly. The index sets for the closed fan Fig. EJa) 

are 



The first row gives the index labels, and the next three rows give the values of the 
indices in the sets /, J, and Jq respectively. Different letters are understood to 
correspond to different values. Additional dots and lines may be present in the 
graph, as long as none of the lines end up in the R-dots shown. These additional 
lines and dots are not drawn because they do not affect the fan relation in any 
way. Similarly, they are not shown in the index sets in (|20|l other than the 



the type [A a ] a [A b ] 




(20) 
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ellipses in the first two rows, which remind us that there may be more lines 
connected to the L-dots. Such ellipses are absent in the last two rows because 
no additional lines are allowed to be connected to the R-dots shown. 

Similarly, the index sets for the fully opened fan, Fig.^b) with all mi = 1, 

are 

... m ... \ 

::: L ) < 21 » 

' ' ' Cm / 

Using S m to denote the symmetric group for the permutation of the first 
m labels, the symmetric groups for Fig. |2Ja) can be read off from (|2U|) to be 
Gi D S m and Gj = Gj Q = S m . Similarly, the symmetric groups for Fig. Hfb) 
can be read off from (|21|l to be Gi 3 S m and Gj = Gj Q = e. We may choose 
Q = e in both cases. Then for Fig. |2fa), TGi — Gi for every T 6 Gj Q , hence 
N[c] = ml J^ReG! ^(-^ e c )- But the ^ as ^ sum ^ s simply the N[c] for Fig. |2fb) 
and l|21|) . Hence it follows from (|16|) that the fan relation (with all mi = 1) is 
true. 

Next, let us derive the double-fan relation l|12(l and i|13|) for the case n a = 
mi, = and m a = rit,. The solution of is now unique. It gives a a = m a = 
rib = rn, and a& = /3 a — Pb — 0. The double-fan relation l|12|) then becomes 

[mAa] =ml[A a ] m . (22) 

The closed double-fan is shown in Fig.[S{a). Its index sets are 




I \ ■ ■ ■ m ■ ■ ■ n ... n + m 
■ b ••• b 



a ■ ■ ■ a 
c ■ ■ ■ c 



(23) 



\ c • • • c 

The opened double- fan is shown in Fig.JSfb). Its index sets are 



label \ 






m 


n 


n + m 






a 


a 


■ b ■ 


b 


J 1 


- 










Jq J 




C ' 




Cl • 


Cm 



(24) 



We may choose Q — (l,n)(2, n + 1) • ■ • (m, n + m) in both cases. For (|23(l . 
Gi 3 S m <8> S' m , Gj — S m , and Gj q = S' m , where S m is the permutation group 
of the first m labels, and S' m is the permutation group for the labels (n, n + 
1, • • • ,n + m). For (PI, Gi d S m ® S' m , but Gj = Gj Q = e. 

For Fig-Bta), TQi = Gi for every TeGj Q . Hence N[c] = ml Y^TeQ, S (Q R e 
c). But the last sum is simply the N[c] of Fig. EJb). In this way 122fl is proven 
by the GTM. 

The fan integral (|10|l can also be obtained from the GTM. It is given by 
Fig. Ufa) without extra dots and lines, or l|2U|) without the ellipses at the end. 
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(a) 



(b) 



Figure 5: Double- fan diagrams with n a 
closed fan; (b) a fully opened fan. 



m& = 0, and m a 



n b 



A 



Then Qj = S m , and N[c] = to! £ fl 5(R 6 c). From (Q2J) and (BJ, we get 

(IJ\IJ Q ) = F(m)=m! ^^—^-^(i?) (25) 



^E? E ^/(^H( fi ) = r (M) 

m!(n — 1)! 



2 



(m + n — 1)! 



(27) 



In getting from i12">l) to i|26[l . the character X* m )(R) = 1 of the totally sym- 
metric representation (to) of the permutation group has been inserted, and 
the orthogonality relation of the characters has been used. To get to (|27f) . 
d(m) = X(m)( e ) = 1 as well as di m -\ = (to + n — l)\/(n — 1)!to! (see l(T5)l ') have 
been used. 

The result in (|27|l agrees with the result (|10|l . It is one of the very few cases 
where the value of the integrals can be obtained group-theoretically without 
knowing the values of the individual £[c]'s. 



3.3 Primitive Integrals 

Integrals in which both symmetry groups Qi and Qj consist only of the identity 
e will be called primitive. This happens when all the indices i a in the set I 
assume distinct values, and all the indices jb in the set J are also different. 
The corresponding diagrams have p dots each on both columns, and precisely 
one solid and one dotted lines connecting to each of the dots. The primitive 
diagrams for p < 3 are shown in Fig. and the ones for p = 4, 5 are contained 
in Appendix IO 

Since Qj = Qj = e, it follows from l|17|) that N[c] = 8{Q € c), where Q 
can be any element of S p . The primitive integrals (|16|l are simply £[c], one for 
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(a) p = 1 



• 

(b) P = 2 



• • * * «t 7^ S» 

• • 4 ^» 

£[(12)(3)] ^[(123)] 
(c) p = 5 

Figure 6: Primitive diagrams for (a) p = 1, (b) p = 2, and (c) p = 3. The 
identity element is everywhere denoted by e. 



each class c of S p . We may therefore use an element of each cycle structure to 
label the primitive integrals, as is done in Fig. |SJ Diagrammatically, the cycle 
structure is translated into the loop structure of its diagram, as can be seen in 
Fig. E| Using (jTHjl along with (fTK|l and the character tables found in Appendix 151 
(note that dj = X/( e ))> the primitive integrals forp < 3 can be easily computed, 
and the results are displayed in Table ^ The results for p = 4, 5 can also be 
found in Table H of Appendix 





c 


p=l 


p = 2 


p = 3 




I 


l 


n 2 -2 


e 


n 


nP-l 


n(n' A — l)(ri 2 — 4) 


(12) 




-1 


-1 




n(n' 2 -l) 


(n a -l)(n 2 -4) 


(123) 






2 






n{n 2 -l){n 2 -A) 



Table 1: Algebraic expressions for the primitive diagrams of p — 1, 2, 3. 
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Figure 7: Arbitrary stack diagram, S(pi,p2 ; ■ ■ ■ ,Pt), of degree p = Yll=i Pi 



3.4 Stack Integrals 

The stack diagrams (see Fig. are direct integrals made up of disconnected 
lines of arbitrary multiplicities. As such, Q = e, and J differs from / only 
by relabelling. Using item 3 of Sec. 2.2, we may assume J = I. Hence stack 
integrals are integrals of the form (II\II). 

Let pi,p2, . . . ,pt be the multiplicities of the disconnected lines in a stack 
diagram. Then Qi = Qj = Q = S Pl ® S P2 &> • • • <8> S Pt , and N[c] is nonzero only 
when the class c is a direct product of the classes Cj of the groups S P( . In that 
case, 

t 

N[c\=Y[ Pl \n t { Cl ), (28) 

i=l 

where ni(a) is the number of elements of S Pi in the class Ci. In other words, 

»*(«) = li: , )f! - { - !)} 

where the class consists of cycles of length j. Denoting the stack integral 
by E(pi,p 2 , . . . ,p t ), we get 

3(pi,P2, ■■-,Pt) = X! n ^ Ci ) ) ® C2 ® " ' ® c *)' ( 30 ) 

Cl,C 2) ... \i=l / 

All stack diagrams can be obtained by making the assignment fi — > pi from 
each representation. In this way, we expect a same number of stack digrams 
as of primitive diagrams, or classes. Using the £ expressions obtained in the 
preceding subsection, the stack integrals for p < 3 can be computed to yield the 
expressions in Tabic [21 
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S(pi, . . . ,p p ) 



sri) = i 

v ' n 



~~ n(n+l) 

2(1,1) = ^ 



5(3) 



3! 



n(ra+l)(ra+2) 



"( 2 ; 1) — (n-l)n(n+2) 
a (l> 1; 1) = „(„2™i)(^2_4) 



Table 2: Algebraic expressions for the stack diagrams of p = 1, 2, and 3. 



3.5 Special Double- Fan Integrals 

The index sets for the fully opened double- fan integrals LA a ] a [Ab] 
with N = 2a) are 

/ label's. 
J 

V J Q I 



1 


a 


a+l . 


. 2a 


b 


■ b 


a 


a 


Cl 


c a 


Ca+l 


■ c 2q 


Q + l 


■ c 2a 


Cl 


c a 



(Fig.Hb) 



(31) 



Hence both Qj and Qj Q consist only of the identity e. As for Qi, it is given 
by S a <E> S a , where the permutation groups S a act respectively on the b and a 
indices in /. The element Q maps Jq to J, i.e. Q = (1, a + l)(2, a + 2) • ■ • (a, 2a). 

The fully opened integral can be computed using (|16fl . with iV[c] given by 
lfT7 | or P ^ . Thus, 



ReQi Tee., 



]T 5(Qi? ec) = X>(3' ec )> 

-ReSj Q' 



(32) 



where the last sum is over every permutation Q' that sends all b indices in 
H31(l to the positions labelled from a + 1 to 2a, and similarly a indices to 
the positions labelled from 1 to a. As a consequence, the allowed cycles of Q' 
must be of even length, and they can be specified by a sequence of nonnegative 
integers (k) = (k\k2 ■ ■ ■ k a ), ki being the number of cycles of length 2i. The 
number of Q' with the class structure (k) that is related to c is given by 



Mel 



(a!) 



n2 



n; 



fa! 



(33) 



In order to see how this is arrived at, consider an example where k\ = 2, k 2 = 2, 
and all other ki values are zero. Then Q' is of the form (ba)(ba)(baba)(baba), 
where the b and a letters should take the distinct index labels in (1, . . . , a) and 
(a + l, ... , 2a) respectively. Another Q' with the same cycle structure can thus 
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a 


[A a ] a [A b ] a 


1 

2 
3 


-l 

n(n 2 — 1) 
2 

(n 2 -l)n 2 (n+2)(n+3) 
-6 


(n-l)n 2 (n+l) 2 (ri+2)(n+3)(n+4)(n+5) 



Table 3: Values of the monomial integrals [v4 a ] Q [^4&] Q for a = 1,2, 3. 




Figure 8: Non-orderly integrals of p = 3. 



be obtained by permuting individually all the a and b labels. This accounts 
for the numerator in (|3~3"jl . However, such permutations do not necessarily give 
distinct Q' elements. The cyclic nature of a cycle tells us that each cycle of 
length 2i will appear i times; this accouts for the i ki factor in the denominator. 
Moreover, no new Q' is obtained if we permute cycles of the same length; that 
accounts for the other factor fcj in the denominator. 

We may now return to (|16fl to calculate the integral [A a ] a [A(,] Q in terms of 
the primitive integrals £[c]. The result for the first few a values are listed in 
Table El 



3.6 Non-Orderly Integrals 

All integrals with degree p < 3 are orderly. The non-orderly integrals of p = 3 
are shown in Fig. |H1 and those related to them by the fan relation ijHJ). The 
calculation of N[c] and the integral for each of them is discussed below. The 
integrals will be labelled by their figure, e.g., integral /(JSJi). 
The index sets for Fig. IBfa) are 



/label\ /l 2 3 

/ 6 b a 

J I I d c c 

\ Jq / \d c c 

They give rise to the symmetry groups Qi = {e, (12)} and Qj — Qj Q — {e, (23)}. 
Moreover, the element Q can be taken to be the identity element. In order to 
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obtain the coefficients N[c] of equation l|16|). we need to compute QTR for all 
T e g jQ and i? e Gi- That triple product is QQ Jq Gi = {e, (12), (23), (132)}. 
As a result, 



J©0 = Z(l, 1, 1) = m + 2£[(12)(3)] + ^[(123)] = 
In the same way, the index sets of Fig. |S[b) are 



(n 2 - l)(n + 2) 



(35) 



(36) 



and hence the symmetry groups are Qj Q = {e} and Gi — {e, (23)}, and the 
exchange element is Q = (12). We thus obtain QGj q Gi — {(12), (123)}, from 
which 





/©)=£[(12)(3)]+£[(123)] = 

follows. 

Finally, for Fig.|^c), the index sets are 



1 



(n 2 - l)n{n + 2) 



(37) 





(38) 



The relevant symmetry groups are Gj Q = {e, (13)} and Gi — {e, (23)}. With 
Q = (12), the set QGj q Gi is {(12), (13), (123), (132)}, and formula (Q2J gives: 

-2 



/©0 = 2(£[(12)(3)]+£[(123)]) 



(n 2 - l)n(n + 2) 



(39) 



The calculation of QQj q Gi is not that cumbersome for p = 3, but it gets 
worse pretty quickly as p increases. For example, let us look at some examples 
of p = 4. 

Let us first calculate Z(2, 1, 1) of Fig. ^b), whose index sets are 



(40) 



Then Q = e, G,j Q = {e, (12), (34), (12)(34)} and Gi = {e, (12), (13), (23), (123), 
(132)}. Thus QGj q Gi = {e, (12), (13), (23), (123), (132), (12), e, (132), (123), 
(23), (13), (34), (12)(34), (143), (243), (1243), (1432), (12) (34), (34), (1432), 
(1243), (243), (143)}, hence 

Z(2, 1, 1) = 2f[e] + 8f[(12)(3)(4)] + 8f[(123)(4)] + 2f[(12)(34)] + 4^(1234)] 

" (41) 





(n- l)n(n + 2)(n + 3) 
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Figure 9: Sigma integral with all lines being of multiplicity one. 



Our last example is the S-integral, shown in Fig. |5J Its index sets are 

/l 2 3 4' 

b b a a 

e d d c 

V e d d c, 




(42) 



Hence Q = e, Gj Q = {e, (23)} and & = {e, (12), (34), (12)(34)}. Multiplying 
such elements accordingly, the set {e, (12), (34), (12) (34), (23), (132), (234), 
(1342)} is obtained for QGj q Qi- Hence 



E = e[e]+3C[(12)(3)(4)] + 2f[(123)(4)]+e[(12)(34)]+C[(1234)] 

71+1 



(n- l)n 2 (n + 2)(n + 3) 



(43) 



4 Comparison of the IM and the GTM 

We have discussed the computation of U(n) integrals (0) in two ways: the IM in 
Sec. 2, and the GTM in Sec. 3. Each of these two methods has its own merits, 
and drawbacks, and in a way they complement each other. The purpose of this 
section is to compare their relative stong and weak points. 

The IM is based solely on the unitarity condition and the invariance 
of the Haar Measure J3J). The method is simple because there is no need to 
know group theory. The conditions relate integrals of the same degree, and 
also integrals of degree p to integrals of degree p — 1. Through these relations, 
specific integrals such as the fan integrals 110(1 and the ^-integrals 10 can be 
obtained, and general relation such as the single-fan relation (jHJ and the double- 
fan relation (|12|l can be worked out. 

The GTM has the advantage of being general, in the sense that all integrals 
can be computed using the formula l|14|) or l(16|) . The price to pay is that we 
have to know the characters of the irreducible representations of the appropriate 
symmetric group, and a triple sum has to be carried out, which can prove to be 
very tedious for integrals of high degrees. Furthermore, unlike the IM, relations 
between integrals are hard to come by, so one must calculate the integrals one 
by one. There are however certain class of integrals, the orderly integrals, for 
which relations can be developed, and the quantity N[c] in (|16fl can be relatively 
easily computed. Then we merely have to know the primitive integrals £[c] in 
1(16(1 to get the value of the orderly integral on hand. The stack integrals 1(30(1 
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(n-3) 



+ 2 




Figure 10: Unitarity sum relation involving the E diagram (second from the left). 
The sum is performed on the index indicated by an arrow. Using the fan relation 
©, the unitary sum can be written as: (2=2 + \)Z{2, 0, 2) + 2E = \Z% 0, 1). 



and the special opened double-fan integrals are examples of this kind. The 
primitive integrals £[c] themselves must be calculated using lJT5)l. 

To summarize, the IM gives a huge number of relations but it is not easy 
to obtain the value of any specific integral. The GTM allows us to calculate 
any specific integral, albeit rather tedious at times, but it is difficult to ob- 
tain relations between integrals. In the next section, we shall discuss a hybrid 
method which makes use of the advantages of both methods. We shall use the 
general GTM formula to calculate a specific set of integrals, and then use the 
IM relations to obtain all the other integrals. 

In the rest of this section, we shall enlarge these general remarks about the 
IM and GTM, by using specific examples presented in the last two sections as 
concrete illustrations. 

The single-fan relation JSJ can be obtained by both the IM and the GTM. 
However the double-fan relation l|12fl in its general form can be obtained only 
by the IM, because most of the integrals involved are not orderly, making it 
hard to derive relations using the GTM. Nevertheless, in special cases involving 
only orderly integrals, l|22J) . the GTM can also be used to derive the relation. 

The Z-formula @ is obtained using the IM, by a series of relations connect- 
ing it down to J dU = 1. Since the Z-integrals are not orderly, it is hard to 
compute them using the GTM except at low degrees. The calculation of those 
by the GTM is shown in equations (|35Jl and 1)41(1 . 

However, since the values of the integrals in the IM are obtained only through 
relations, it may be relatively complicated to calculate just one specific integral. 
This is where the GTM is superior, because of the general formula (|14|) valid for 
any one integral. For example, it is easy to obtain the E integral (|43l) . assuming 
of course the £[c]'s to be already known. We can also obtain it using the IM, 
as we shall show below, but that involves a few steps because we must get it 
from relations. To see how that is done, look at Fig. which is the unitarity 
relation applied to the dot of the first diagram indicated by an arrow. The first 
and third diagrams can be related to Z(2,0, 2) by using the fan relation JSJ), 
and similarly the diagram on the right can be related to Z(2, 0, 1). Using the 
Z-formula 0, we then obtain 

- - \ 



2(2,0,1)- (j^ + A ^(2;0,2) 
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n + 1 

(n- l)n 2 (n + 2)(n + 3)' 
the same as the result (g3Jl obtained by the GTM 

5 Hybrid Method 

Having understood the relative merits of the GTM and the IM, it is possible to 
combine their strengths into a more efficient hybrid calculational scheme. 

The strategy is to start with one or more integrals that can be computed 
by the GTM with relative ease. Generally speaking, such integrals are ordered. 
Once they are obtained, the many relations of the IM can be used to calculate 
other integrals from them. 

To illustrate this strategy, we will consider how the hybrid method can be 
used to calculate all double- fan integrals. 

By a double-fan integral, we mean any integral with two L-dots and any 
number of R-dots. Fig. |3[a) shows a closed (double-fan) integral (with the 
understanding that there are no extra dots or lines than those shown), and 
Fig- Bib) shows a fully opened (double-fan) integral. We may also have partially 
opened (double-fan) integrals, in which every branch, namely, every R-dot with 
its connecting lines, can be regarded as a closed integral. See Fig. for an 
example of a partially opened integral. 

As in Sec. 2.2, a fully opened integral is denoted by [A a } a ° [A b } a " [B a f° [B b f\ 
and its corresponding closed integral is denoted by [a a A a + a b A b + f] a B a + (3 b B b ]. 
For a partially opened integral, we will denote it as a product of the closed 
integrals of each branch. See Fig. ^] for examples. 

Using Ijl2|l and all double- fan integrals can be expressed as sums of 

fully opened integrals. Integrals of the form [A a ] Q [yl( ) ] Q are given by and 
Table|31 but we still have to know how to calculate a fully opened integral when 
Pi ^ 0. As shown in Appendix [D] the IM allows us to relate them to those with 
fti = 0, by using the following formula 

[A a ] a [A b ] a [B a f"[B b f<> = 

min(/3 a ,/3f,) , , \ / \ 

£ j(-ire!(^j (^ b j(n + 2 a -l + 2e). (45) 

(n + 2a-2 + e)l(n + 2a-l + 2e)\ + +e l 

(n + 2a + l3 a -l + e)\{n + 2a + l3 b -l + e)\ [ aJ [bJ J ' 

We close this section by showing how to use (|12|l and to calculate the 
integrals in Fig. 

5.1 Fig. [TT](a) 

There are two equivalent forms for this diagram. One is 

[A a + 2A b ][A a ] = (2[A a ][A b } 2 ) [A a ] = 2[A a ] 2 [A b } 2 (46) 
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(44) 




(a) (b) 

Figure 11: Partially opened double- fan integrals, (a) There are two equivalent 
forms for this graph: [A a + 2Ab] [A a ]&nd [Ab + B a + Bb] [A a ], (b) There are four 
equivalent forms for this graph: [A a + Ab + B a ] [A a + Ab] , [2B a + Bb] [A a + Ab] , 
[A a + A b + B a ] [B a + B b ] , [2B a + B b ] [B a + B b ] . 



and the other is 

[A b + B a + B b ][A a ] = 4 ([A b ][B a ][B b ]) [A a ] = A[A a ][A b ][B a ][B b ], (47) 

where equations (|12|) and l|13l) have been used. The integral I ljllb ) is obtained 
by adding up gBJ and gZJ. 

Using l|45l) . we can express all fully opened integrals in the form [A a ] Q [Af,] a . 
Applying to the present case, we get 

[A a ][A b ][B a ][B b ] = —L^[A a ][A b ] -^—[A a f[A b ] 2 . (48) 
[n + l) (n + I) 

Using Table |3 we finally obtain 

JCQb) = 2[A a ] 2 [A b ] 2 + A[A a ][A b ][B a ][B b ] 



■In r . , 2 u v> I 

-4 

(n 2 - l)n(n + 2)(n + 3)' 



5.2 Fig. [TT](b) 

As shown in Fig. Illf b). Jljllt)) has four equivalent forms. For one branch, the 
factors are 



[A a + Ab + B a ] = 4[A a ][A b ][B a ], 

[2B a + B b ] = 2[B a ] 2 [B b ]; (49) 

and for the other branch, they are 

[A a + A b ] = [A a ][A b ], 

[B a + B b ] = [B a ][B b ]. (50) 
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Hence 



JdW = 4 [^a] 2 [A b ] 2 [S a ] + 6 [A a ][A b ][B a ] 2 [B b ] + 2 [B a f[B b ] 2 . (51) 

We will now express each of the three monomial integrals in l|51|l in terms 
of L4 a ] a L4 b ] Q . First, with respect to g3J), L4 Q ] 2 L4 b ] 2 [_B a ] is characterized by 
a = 2, j3 a = 1, and j3 b = 0. The vanishing of /? b causes ((45(1 to consist of the 
single term: 

[A a ] 2 [A b ] 2 [B a ] = -^[A a f[A b f. (52) 
(n + 4) 

Second, L4 a ][^l&][-Ba] 2 [-B&] has a = 1, a = 2, /3& = 1, and the sum in l|4*5*|l gives: 



[A Q ][A b ][£ Q ] 2 [i? b ] = - ; — ■ ; — [Aa] [Ab] - ; ^ ; ^ L4 a ] 2 [A b ] 2 . (53) 



(;; -h2) 2 (n + 3)^°^ (n + 2)(n 4) 
Finally, [S a ] 3 [_B b ] 2 , having a = 0, (3 a = 3, /3 b = 2, can be expressed as 



[Ba]W _ 7i 2 (n + l) 2 (n + 2) n(n + 2) 62 (n + 3) [AanA] + 

° [A a ?[M 2 (54) 



(ra + l)(n + 2)(n + 4) 



from equation 1)45(1. Using the fan relation, notice that [i? a ] 3 [i? b ] 2 can also be 
reduced to ^Z(3,0,2). 

The expressions of [j4 a ]L4 b ] and L4 a ] 2 L4 b ] 2 in terms of n have already been 
determined in Example 1. The final answer is obtained by inserting l|52 |l -l|54 |l 
into 1)51(1 . The result is: 

4( " 2+2) MA* 



(n + l)(n + 2)(n + A) 
2(n 2 + 2n + 4) 



(n 2 - l)n 2 (n + 2){n + 3)(n + 4) ' 
which can be verified using the plain group theoretical formula 1(14(1 . 



6 Conclusion 

In this article, we have pursued the goal of finding an efficient method to calcu- 
late the monomial integral or J5J). We find that the IM discussed in Sec. 2 
is superior for deriving relations between integrals, but the GTM is able to give 
a formula to calculate any integral. The GTM formula involves a triple sum 
whose computation is often tedious and prone to mistakes. The sums simplify 
for orderly integrals, in which the invariant groups Qi and Qj Q are either dis- 
joint, or one is contained in the other. For non-orderly integrals, the hybrid 



21 




(c) (d) 

Figure 12: Diagrams that result from rotating away the pairs of lines (a) [A a ], 
(b) [Af,] , (c) [B a ] , and (d) [Bf,] , from the compound pattern of Fig. [3J 

method is probably the most efficient. It uses the IM to relate them to some 
orderly integrals that can be calculated by the GTM with relative ease. 

This research is supported by the Natural Sciences and Engineering Research 
Council of Canada and by the Fonds de recherche sur la nature et les technologies 
of Quebec. 



Appendices 
A Derivation of Equations (1121 ) and (II 3D 

To prove and l|13fl . we use the rotation technique discussed in item [5] of 
Section 12.11 and equation Q , to spin off from the R-dot of Fig. IHta) a new 
R-dot attached to a pair of solid-dotted lines. Depending on whether the basic 
pattern of this new R-dot is [A a ], [At], [B a ], or [_£?&], we get the graphs shown 
in Fig. Eta), Efb), Etc), andEJd) respectively. 

Repeating this spin-off operation over and over again on the R-dot still 
containing a compound pattern, eventually we come to a graph where every 
R-dot is given by a basic pattern. The resulting integrals are given in (fL2"|) . 
corresponding to the decompositions of the compound pattern into a sum of all 
possible fully opened integrals obtained by spinning off. 
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The coefficient v(a' a a' b f3' a f3' b ) of these basic integrals is derived from a com- 
bination of three factors: 

1. Each time that we spin off a basic pattern from a compound pattern with 
d pairs of solid-dotted lines, there is a factor 1/d arising from eq. J7J), by 
taking e = 1. Since we start from a compound pattern with N pair of 
lines, by the time we come to a fully-opened pattern we have accumulated 
a factor 1/AM. 

2. The TV R-dots in the final pattern that is fully opened can be spinned off 
in a different sequential order. According to J7J), they must be summed 
over. This gives rise to a factor N\/a' a \a' b \{3' a \[3 b \. 

3. At any time when we spin off a basic pattern, we can choose its single pair 
of solid-dotted lines in all possible ways. Eq. Q says that we must sum 
over all these possibilities. The mulitipliciy factor is given by the number 
of permutations of these lines that lead back to the same basic integral. 
It is a factor of m a \ mj,! nj n^. 

Assembling these three factors, we get 

v(a a a b( 3 M = a , !a , !/w , (55) 
which is identical to (|13[l . 

B Character Tables 

The character tables for p — 2,3 are given here in the form used by M. Hamer- 
mesh in fSJ. The rows are labelled by the partitions that define the represen- 
tations, and the columns are labelled by the cycle structures that define the 
classes. The number of elements in each class, n(c), is written above the classes. 
The table for p = 1 is trivial, and it consists of the sole value 1. 

B.l p = 2 



Part. \ Class 


1 1 

(I 2 ) (2) 


(2) 
(I 2 ) 


1 1 

1 -1 



B.2 p = 3 





1 


3 


2 


Part.\Class 


(I 3 ) 


(1,2) 


(3) 


(3) 


1 


1 


1 


(2,1) 


2 





-1 


(I 3 ) 


1 


-1 


1 
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C Primitive Diagrams for p = 4 and p = 5 



Using l|14f> and the character tables for S4 and S5 [5j , the algebraic expressions 
for the primitive diagrams of Fig. 1131 and Fig. 1141 can be obtained, and they are 
given in Table 0] 



£H £[(12)(3)(4)] £;rU23X4)l £[(12)(34)1 £[(1234)1 

Figure 13: The p — 4 primitive diagrams. 



£[*] 



£[(12)(3)(4)(5)] £[(123)(4)(5)] £[(12)(34)(5)] 



£[(1234)(5)] £[(123)(45)] £[(12345)] 

Figure 14: The p — 5 primitive diagrams. 
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?(cq) 




p — 4 


p = 5 




n 4 -8n 2 +6 


n 4 -20n 2 + 78 


e 


n 2 (n 2 -l)(n 2 -4)(n 2 -9) 


n(n 2 - 


-l)(n 2 -4)(n 2 -9)(n 2 - 


16) 


(12)(3)(4) 


-1 




-(n 2 -2)(n 2 -12) 




n(n 2 -l)(?i 2 -9) 


n 2 (n 2 


-l)(n 2 -4)(n 2 -9)(n 2 - 


-16) 


(123)(4) 


2n 2 -3 




2 




n 2 (n 2 -l)(n 2 -4)(n 2 -9) 


n( 


n 2 -l)(n 2 -4)(n 2 -16) 




(12)(34) 


n 2 +6 




n 2 -2 




n 2 (n 2 -l)(n 2 -4)(n 2 -9) 


n(n 2 - 


-l)(n 2 -4)(n 2 -9)(n 2 - 


16) 


(1234) 


-5 




-5n 2 +24 




n(n 2 -l)(n 2 -4)(n 2 -9) 


n 2 (n 2 


-l)(n 2 -4)(n 2 -9)(n 2 - 


-16) 


(123)(45) 






-2(ri 2 +12) 






n 2 (n 2 


-l)(n 2 -4)(n 2 -9)(n 2 - 


-16) 


(12345) 
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n(n 2 - 


-l)(n 2 -4)(n 2 -9)(n 2 - 


16) 



Table 4: Algebraic expressions for the primitive diagrams of fourth and fifth 
degrees. In the p — 5 case, the elements from row two to row five should be 
written with the additional (5) one-cycle. 



D Derivation of Equation ( 1451) 

We would like to express the general fully-opened integral [^4 a ] a L4;,] a [i3 a ] /3a [-Bb]' 3 '' 
in terms of the special ones of the form [A a ] a [A b ] a . The idea is to apply a uni- 
tarity sum on the [B a ] or [B b ] basic patterns to get rid of them. To get the 
final result we also need to apply the fan relation (JHJ) or the double-fan relation 
l|12|) and (|13|) . Our approach is to first determine how can [A a ] a [Ab] QI [B a ]^ a 
be reduced to fully opened integrals involving only the [A a ] and [Ab] patterns. 
With such an information at hand, we will then try to reduce the more general 
[A a } a [A b ] a [B a f"[B b f» integrals into the [A a } a [A b ] a [B a f° integrals. 

Let us apply a unitarity sum on one of the [B a ] patterns in L4 a ] Q L4;,] Q [i? a ] /3a : 

(n - (2a + (3 a - 1)) [A a ] a [A b ] a [B a ]^ + [fi a - 1) [A a ] a [A b ] a [B a ]^- 2 [2B a ] 
+ a ( [AXlAbr^iAi + B^lBaf'- 1 + [ij " 1 ^ + BJ^]"^]"-- 1 ) 
= [A a ] a [A b ] a [B a f°- 1 (56) 

Using (jHJ), [2£? a ], in the second term above, can be rewritten as 2[£? a ] 2 . Fur- 
thermore, equation l(T3|l tells us that [A b + B a ] and [A a + B a ] can respectively 
be rewritten as 2L4f,][i? a ] and 2L4 Q ][£? a ], and the term in parentheses above 
simplifies to 4L4 a ] Q L4fc] Q [£? a ] / V As a result, relation (|56|l reduces to: 

[A a r[A h r[B a f° = 1 — - [A a rw{B a }^-\ ^ 

n + 2a + p a — 1 

Using the relation l|57|l recursively on its right-side, until no [B a ] remains, we 
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obtain: 



[A a riA b } a {B a f« 

(n + 2a + J a - 1) (n + 2a + J a - 2) " ' (n + 2a) 
(n + 2a- 1)! 



{n + 2a + (3 a -l)\ 



[A a ] a [A b ] a , (58) 



and the first step of the work is completed. 

Assuming that (3 b < /3 Q , let us perform a unitarity sum on a pattern in 

[A a ] a [A b } a [B a }^[B b ]^: 

(n-(2a + (3 a + p b -l)) [A a ] a [A h ] a [B a f° [B b f> 
+ (J3 b - 1) [A a } a [A b } a [B a }^[B b }^[2B b } 
+ (3 a [AXWlB^-^Ba + B b \[B b \ 0b ~ 1 
+ a ( [A a ] a {A,}*- 1 [A b + B b ] [B a ]P> [B b f^ 

+ [AX' 1 [A a + B b ] [A b ] a [B a f* [B^- 1 ) 
= [AX[A h r[B a f"[B h f»-\ (59) 

Relation ©, or formula l|13|> . again permit to make some simplifications, i.e. 
[2B b ] = 2[B b } 2 , [B a +B b ] = [A a ][A b } + [B a }[B b ], [A b +B b ] = 2[A b ][B b ], [A a +B b ] = 
2[A a ][B b ]. By making the proper substitutions in l|59(l . the recursion equation, 

[A a ] a [A b ] a [B a f°[B b f> = 

{[A a ] a [A b ] a [B a f«[B b ]^- 1 



{n + 2a + (3 b - 1) 

- (3 a [Aar+^Abr+^Baf^iBb]^- 1 } , 



results. The above can be solved to give 

[A a } a [A b ] a [B a }^[B b f b = 

'M {ft 



Hb f 



(n + 2a - 1 + 2e) 



(n + 2a-2 + e)! r ^ + + A _ e 



(n + 2a + /3 h - 1 + e)! 
which upon substitution of yields the desired equation (|43)l 
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Abstract 

In a previous article, an 'invariant method' to calculate monomial 
integrals over the U(n) group was introduced. In this paper, we study the 
more traditional group-theoretical method, and compare its strengths and 
weaknesses with those of the invariant method. As a result, we are able to 
introduce a 'hybrid method' which combines the respective strengths of 
the other two methods. There are many examples in the paper illustrating 
how each of these methods works. 

1 Introduction 

This article deals with the calculation of integrals of the form 

J(dU)U* ljl ---U* pj U klh ---U kqlq (1) 

over the U(n) group, where (dU) is the invariant Haar measure normalized 
to J(dU) = 1, and Uij is a U(n) matrix element, with U*j being its complex 
conjugate. 

These integrals and their generating functions are useful in many areas of 
physics, including two-dimensional quantum gravity [1], QCD, matrix models, 
and statistical and condensed-matter problems of various sorts [2]. They are 
also needed in the parton saturation problem at small Feynman-x [3] . 

The integral (1) depends on the indices I = - ■ ■ i p }, J — {ji ■ ■ -j p }, K = 
{k\ ■ ■ -k q }, and L = {l\ ■ ■ -l q }, so it will be denoted as (IJ\KL): 

(IJ\KL) = J (dU)U!jU KL , (2) 

where Ufj — Y\ P a= \U* a j a , and similarly for Ukl- Since the matrix elements 
commute, Ufj — Uf pJp where Ip — {ip(i) • • - ip( p )} is obtained from / by a 
permutation P £ S p of its p indices. Hence 

(IJ\KL) = (I P Jp\K T L T ) (3) 
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for any P € S p and T e S q . 

The integral is nonzero only when p = q, a number which will be referred 
to as the degree of the integral. Without loss of generality, it turns out that we 
may assume K = I and L to be a permutation of J, namely, L = Jq for some 
Q E S p . The value of the integral depends on what the index sets /, J, and 
what the element Q are, so even for a given p, there are many distinct cases. 
The best way to distinguish them is to represent each integral by a diagram in 
a way to be explained in the next section. 

Integral (1) has been computed using a graphical technique [4]. It can also 
be obtained using the Itzykson-Zuber formula [5] as a generating function, or 
directly from group theory [6] using the Frobcnius formula [7] . We shall refer 
to this last method as the group-theoretical method, or GTM for short. In the 
GTM, a general formula is available to compute (1). It involves a triple sum 
over an expression containing characters of the symmetric group S p , as well as 
the dimensions of irreducible representations of S p and of the unitary group 
U(n). One of the sums is taken over all the relevant irreducible representations, 
and the others are taken over the symmetry groups of the index sets / and J. 
These sums could be long and tedious for a large p, and for most symmetry 
groups. 

A different way to calculate (1) was introduced in a recent paper [8]. This 
method relics only on the unitary nature of the matrix elements, and the in- 
variance of the Haar measure. In particular, no knowledge of group theory is 
necessary. The invariance of the Haar measure, as well as the off-diagonal uni- 
tarity relation, are used to derive relations between integrals of the same degree. 
The diagonal unitarity relation connects integrals of degree p with ones of degree 
p — 1. Through a chain of these relations, the desired integral is finally related 
to the basic integral of degree 0, which is / dU = 1. The desired integral is then 
solved from this chain of relations. We have called this method the invariant 
method, or IM for short. 

The purpose of this paper is to compare the pros and cons of the GTM 
with the IM. In order to do so we must first study and understand better the 
nature of the GTM. Armed with this comparison, we will be able to design a 
new hybrid method which combines the strengths of these two other methods. 

The IM is reviewed in Sec. 2. It is used to derive a new 'double-fan' relation 
needed in the example for the hybrid method. 

The GTM is reviewed and studied in Sec. 3. It can be used to derive simple 
relations, but the relations derived in this way are nowhere as powerful as those 
derived with the IM. The group-theoretical formula can be used to calculate 
any integral, but generally that is tedious and has to be done integral by in- 
tegral. However, for integrals whose two symmetry groups arc disjoint, or one 
is contained in the other, systematics emerge to make the calculation simpler. 
Several of these 'orderly' integrals are studied in Sec. 3. 

A comparison of the strengths and weaknesses of the two methods is to be 
found in Sec. 4. Armed with an understanding of their relative merits, we design 
a 'hybrid method' in Sec. 5 to take advantage of their respective strengths. This 
method is illustrated by the calculation of a class of 'double-fan' integrals. There 
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are also four Appendices showing details of various calculations. 

2 The Invariant Method (IM) 
2.1 A Brief Review 

The invariant method presented in a previous paper [8] can be used to calculate 
the integral (1). The method exploits the unitarity of the U(n) group elements, 

n n 

E^ = E^« = ^ ( 4 ) 

and the invariance of the Haar measure, in the form 

y (d£7) f(U, U*) = j (dU) f(U*, U) = J (dU) f(U T , U* T ) 

= J(dU)f(VU,V*U*) = J(dU)f(UV,U*V*), (5) 

for any function / and any V £ U(n). 

The unitarity relation (4) relates integrals of the same degree if i ^ I, and it 
relates integrals of degree p to integrals of degree p — 1 if i = I. Other relations 
between integrals of the same degree can be obtained from (5), by suitable 
choices of V. Here are some of them discussed in the previous paper [8]. 

1. Using Vij = e^Sij, it follows that p must be equal to q in order to avoid 
the vanishing of integral (2). The number p will be called the degree of 
the integral. 

2. Using Vij — e lc ^ i 5ij, it follows that (2) is nonzero only when K = Im and 
L = Jr for some M,R £ S p . Using (3), and denoting i?M _1 by Q, only 
integrals of the type 

(IJ\IJ Q ) = J(dU)UfjU IjQ (6) 

are nonzero, so from now on we need to consider only integrals of this 
type. 

An integral with Jq — J will be called a direct integral. For such integrals 
we may always choose Q = e, the identity permutation. Otherwise, the 
integral is an exchange integral. 

Integrals are represented diagrammatically as follows. Each distinct value 
in the index set / is represented by a dot on the left (L-dot), and each 
distinct value of the index set J or Jq is represented by a dot on the right 
(R-dot). The factor U*j is shown as a thin solid line between the L-dot 
i and the R-dot j, and the factor Uij is shown as a dotted line between 



3 




Figure 1: Examples of U(n) integral diagrams, (a) The unique exchange integral 
for p = 2; (b) a Z-integral with arbitrary multiplicities mi, to 2 , and TO3. 



these two dots. The factor U*^Uij is represented by a thick line, or more 
generally, the factor U*™U™j is represented by a thick line with a pair of 
numbers (to, n) written beside it. If to = n, then only a single number to 
is written. The numbers (to, n) or to will be known as the multiplicities 
of the line. See Fig. 1 for an illustration. 

3. With V chosen to be a permutation matrix of n objects, it follows that 

(IJ\KL) = (I'J\K'L) = {IJ'\KL'), 

where I' is obtained from I by a reassignment of the values of its indices, 
e.g. if / = (334), then I' may be (558). K 1 is obtained from K by the 
same reassignment, and similarly for J' and V . As a result, there is no 
need to know the values of the indices of the L-dots, nor the R-dots. This 
is why the dots in Fig. 1 are not labelled. 

4. As a consequence of the first two equalities in (5), an integral remains the 
same under the interchange of the solid lines with the dotted lines, or the 
L-dots with the R-dots. 

5. Using V = R(ab), the rotation matrix in the (a, b) plane, a 'spin-off re- 
lation' is obtained. Consider a R-dot imbedded in an arbitrary integral 
M = (IJ\IJq), with d pairs of solid-dotted lines attached to the dot. Now 
spin off e pairs of these lines to create a new R-dot and a new integral. 
There are many ways to choose the e pair of lines, each possibly corre- 
sponds to a different integral. Let M e be the sum of all these integrals. 
Then the quantities M and M e are related by the spin-off relation 

M e = M ^fj, (7) 
where (^^j is the binomial coefficient. 

The relation is local in that it is independent of the the structure of the 
rest of the diagram. The same relation can also be used to spin off a L-dot. 
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(a) (b) 

Figure 2: The fan diagrams shown here can be a part of a larger diagram. In 
that case, there may be many more dots and many more lines in the complete 
diagram, provided none of the additional lines land on the R-dots shown, (a) A 
closed fan; (b) a partially opened fan. If all m, = 1, then it is said to be a fully 
opened fan, or just an opened fan. 



In what follows, we summarize two general results obtained in the previous 
paper [8] using the IM. 

2.1.1 The fan relation 

The fan relation 

(11* m ') r 

(dU)A\U aCl \ 2mi ---\U ac f mt = " (dU)A\U ac \ 2m (8) 

relates integrals of the same m — X)*=i TO i: where A is an arbitrary prod- 
uct of matrix elements of U and U* whose column indices are different from 
Ci, C2, . . . , Ct- The column index c on the right could be taken to be one of the 

Cj's. 

Diagrammatically, the integral on the right of (8) is shown in Fig. 2(a), and 
the integral on the left is shown in Fig. 2(b). The additional lines and dots 
corresponding to the factor A are not shown, because they do not affect the 
spin-off relation. We shall refer to Fig. 2(a) as a closed fan, and Fig. 2(b) as a 
partially opened fan. If every mj = 1 , then it will be said to be a fully opened 
fan, or simply an opened fan. 

In particular, a closed fan integral is ml times an opened fan integral. In 
fact, this relation between the two types of fans immediately gives rise to the 
relation (8) between a closed fan and a partially opened fan. To see it, note 
that each branch of a partially opened fan is itself a closed fan. By opening up 
all of them, we get the fully opened fan integral, multiplied by a multiplicity 
factor Ylj m r from all the branches. Thus a closed fan is to!/]^ rrijl times a 
partially opened fan, as given by (8). 
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(a) 



(b) 



Figure 3: Double-fan diagrams, (a) A closed diagram; (b) a fully opened dia- 
gram. 



2.1.2 The Z-integral and the fan integral 

In [8], we also obtained a general formula for integrals of the type shown in 
Fig. 1(b), for arbitrary nonncgative integers mi,m 2 , and m 3 . We call that the 
'Z-formula' because of the shape of the graph. It is 



which is the integral for Fig. 2(a) when there are no additional dots or lines 
around. 

2.2 Double-Fan Relation 

The fan relation can be generalized to a double-fan relation, connecting the 
closed 'double-fan' diagram of Fig. 3(a) with a (fully) opened double-fan diagram 
such as Fig. 3(b). As in Fig. 2, there may be additional dots and lines in the 
integral, but none of them may end up on the R-dots shown. From that relation, 
we can also deduce relations between a closed double-fan and a partially opened 
double-fan, as done in the single-fan case. 

The double-fan relation is considerably more complicated than the single-fan 
relation (8), because there are many more double-fan graphs. Each R-dot of a 
(fully) opened (double-fan) graph such as Fig. 3(b) falls into one of four basic 
patterns: [A a ], [Ab], [B a ], and [Bb], shown in Fig. 4. If the solid and dotted lines 
end up on the same L-dot, the pattern is a [B]; otherwise it is an [A]. The 
subscripts a and b tell us which L-dot the solid line emerges from. 




■ (9) 



(n + mi — 2)! (n + m.3 — 2)! (n + mi + m 2 + m.3 — 1)! 
In the special case where m 2 = m 3 = 0, this becomes 




(10) 
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a 



a • 



b •' 



b • 



lA a ] 



[A b ] 



[B a ] 



[B b ] 



Figure 4: The four basic patterns for the R-dots of a fully opened double-fan 
graph. 



Suppose there are «i number of [Ai] and Pi number of [Bi\ patterns in a 
(fully) opened (double-fan) graph. Then there are m a solid and n a dotted lines 
emerging from the L-dot a, and m b solid and n b dotted lines emerging from the 
L-dot b, where 

m a = a a + f3 a , n a = a b +{3 a , 

m b = oeb + fib, rib = a a +(3 b . (11) 

The total number of R-dots in the opened graph is ./V = m a + m& = n a + rib- 

When these N R-dots are merged together, we get the closed (double-fan) 
graph depicted in Fig. 3(a), which will be denoted by [(m a n a )(m b nb)]- If N = 1, 
this is just one of the four basic patterns discussed before. If N > 1, we will 
call it a compound pattern. 

From (11), we see that if we replace at and by a[ = ai and fi[ — fa — 
with any integral £ which keeps a\ and [3[ nonnegative, then we get the same 
closed graph by collapsing this new opened graph. Conversely, it will be shown 
in Appendix A that the closed graph [(m a n a )(mbnb)] can be spinned off into 
a sum of several opened graphs, one for each (a-,/3-). The double-fan relation 
expressing that quantitatively is 

[(m a n a )(m b n b )} - ^ v(a' a a' b (3' a (3'b) [A a }< [A b ] a '" [B a ] ^ [B b ] # 

= YlKAa+a' b Ab + p' a B a + p b Bb], (12) 

where 

v(a a a b (3 a (3 b ) = a n a , blp , Jp , r (13) 

and the sum is over all solutions (a' a , a' b , (3' a , (3' b ) of (11). 

The double fan becomes a single fan if the L-dot b is not connected, namely, 
if m b = n b = and m a = n a = m. In that case (12) becomes 

[mB a ]=m\ [B a ] m , 

which is just (8) (when all rrn = 1) in another notation. 
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3 Group-Theoretical Method (GTM) 



3.1 A Brief Review 

Using group theory, the integral (6) can be turned into a multiple sum [5, 6]. 
In the notation used in Appendix A of [8], the formula is 

(u\uq)= E E Ett^x*/^), (u) 

where p is the degree of the integral. The symbols Qi and Qj represent the 
symmetry groups of the row and column index sets. More precisely, Qx = 
{P € S P \P{X) = X}, with X being either I or J. The irreducible represen- 
tations of the symmetric and unitary groups are both labelled by a sequence 
/ = (h, h, f P ), with /i > f 2 > ■ ■ ■ > f p > 0. Xf(P) is tnc character of 
P G S p in the irreducible representation with signature /. The dimension of the 
irreducible representation / is given by df — Xf( e ) f° r S p , and by df for U(n). 
A formula for df is: 

i (f f , D(h + (n-l),f 2 + (n-2),...,f n ) 
d{fl,...,f„) = —, rr , (15) 

D{n — 1, n — 2, . . . , 0) 

where D(xi, . . . , x n ) is the Vandermonde determinant given by ni<fe( x i — x k)- 
Since Xf(g) depends only on the class c g that g belongs to, we may write it 
as Xf( c g)- With this notation, eq. (14) can be re- written as 

(7J|/Jq> = E^[cK[c], (16) 

C 

where the sum is taken over all classes c of 5 P , 

N ® ^ E E 5 ( S Q R e c ) 

ReGi seg., 

is the number of elements of the type SQR in the class c, and 

^ c ] = Et4i-X/(c). (18) 

It is not difficult to see that Q is not unique, because Q' = SQ = QT for 
any S £ Qj and any T e Q Jq is another possible Q. It does not matter which 
Q we pick in (17). That equation can also be written as 

N ® = E E 6 ^ TR e c )- ( 19 ) 

ReQi Teg., Q 

It is straightforward but generally very tedious to compute N[c], because 
we need to calculate the product QTR for every T £ Gj q , every R G Qi, and 
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determine what class c the product belongs to. Then we have to count up all 
the products that are in a given class c to get N[c]. However, the task becomes 
considerably more manageable if cither Qj and Q Jq are disjoint, or if one is 
contained in the other. We shall refer to integrals with those properties as 
orderly. Further simplification occurs for direct integrals, because in that case 
Q can always be chosen to be the identity e, so the triple product is reduced to 
a double product TR. 

The calculation of £[c] in (18) is simpler than the calculation of N[c], but 
still we know of no closed form of it valid for every class c and every symmetric 
group S p . The best we can do is to compute them case by case. Results are 
given in Sec. 3.3. Each £[c] is actually an orderly integral with Qj = Qj = e, to 
be referred to as a primitive integral. 

Other integrals can be computed in terms of the primitive integrals, if N[c] 
is known. We shall discuss two orderly integrals for which N[c] can easily be 
obtained. In Sec. 3.4, we discuss the stack integrals, which are direct integrals 
with Qj = Qj. In Sec. 3.5, we discuss the fully opened double-fan integrals of 
the type [A a ] a [A b ] a . 

Relations between orderly integrals may be obtained without knowing the 
explicit values of £[c], if their iV[c]'s are related in a simple way. This is the case 
for the single-fan relation, and the double-fan relation with n a — ra b = and 
m a = n b = m. They will be discussed in Sec. 3.2. However, general double-fan 
integrals are not orderly, so we cannot obtain the general double-fan relation 
by the group-theoretical method, at least not in the present way. We will also 
show that the closed (single-)fan integral can also be computed without explicitly 
knowing what £[c] are. This is one of the very few cases where integrals can be 
obtained group-theoretically without explicitly knowing £[c]. 

That leaves the non-orderly integrals, for which each term of the summand 
in (19) has to be calculated separately to get N[c\. The first non-orderly integral 
occurs in degree p = 3. In Sec. 3.6, we shall show how to calculate some of the 
p = 3 and p = 4 non-orderly integrals. 

3.2 Single-Fan and Simple Double-Fan Relations 

The single-fan integrals are orderly. The index sets for the closed fan Fig. 2(a) 
are 



/ label \ 




■ m 








a 








c 




\J Q ) 




c 


J 



(20) 



The first row gives the index labels, and the next three rows give the values of the 
indices in the sets I, J, and Jq respectively. Different letters are understood to 
correspond to different values. Additional dots and lines may be present in the 
graph, as long as none of the lines end up in the R-dots shown. These additional 
lines and dots are not drawn because they do not affect the fan relation in any 
way. Similarly, they are not shown in the index sets in (20) other than the 
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ellipses in the first two rows, which remind us that there may be more lines 
connected to the L-dots. Such ellipses are absent in the last two rows because 
no additional lines are allowed to be connected to the R-dots shown. 

Similarly, the index sets for the fully opened fan, Fig. 2(b) with all m, = 1, 

are 



/ label \ 
/ 
J 

V J Q J 



Z 1 

a 

Cl 
Cl 



to 
a 

Cm 



J 



(21) 



Using S m to denote the symmetric group for the permutation of the first 
to labels, the symmetric groups for Fig. 2(a) can be read off from (20) to be 
Qi D S m and Qj = Q,j Q — S m - Similarly, the symmetric groups for Fig. 2(b) 
can be read off from (21) to be Qi D S m and Qj = Q Jq = e. We may choose 
Q = e in both cases. Then for Fig. 2(a), TQi = Qi for every T e Qj Q , hence 

N[c] = m! Y^Reg, S ( R e c )- But tnc last sum is sini p!y tnc N [ c \ for Fi S- 2 ( b ) 
and (21). Hence it follows from (16) that the fan relation (with all rrii — 1) is 
true. 

Next, let us derive the double-fan relation (12) and (13) for the case n a = 
rah = and m a = rib. The solution of (11) is now unique. It gives a a = m a = 
rib = m, and ab = (3 a — fib — 0. The double-fan relation (12) then becomes 



[mA a ]=m\ [A a ] m . 
The closed double-fan is shown in Fig. 5(a). Its index sets are 



(22) 



label \ 






• TO 


■ n 


. n + to 






<\ : 


a 


■ b ■ 


b 


>) 


- 




c 






Jq 1 








c ■ 


c 



(23) 



The opened double-fan is shown in Fig. 5(b). Its index sets are 



label \ 






TO 


n 


. n + to 










a 


■ b ■ 


b 
















Jq / 




c ■ 




ci • 




J 



(24) 



We may choose Q = (1, n)(2, n + 1) • • • (to, n + to) in both cases. For (23), 
Qi D S m ® S' m , Qj = S m , and Q Jq = S' m , where S m is the permutation group 
of the first to labels, and S' m is the permutation group for the labels (n, n + 
1, • • • , n + to). For (24), Qj D S m <g> S' m , but Q.j = Q. Jq = e. 

For Fig. 5(a), TQi = Qi for every T e Qj Q . Hence N[c] = to! Etgs, 5 (Q R e 
c). But the last sum is simply the N[c] of Fig. 5(b). In this way (22) is proven 
by the GTM. 

The fan integral (10) can also be obtained from the GTM. It is given by 
Fig. 2(a) without extra dots and lines, or (20) without the ellipses at the end. 
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(a) 



(b) 



Figure 5: Double-fan diagrams with n a = m b = 0, and m a = n b = m. (a) A 
closed fan; (b) a fully opened fan. 



Then Gi = S m , and N[c] = m\ Y,r S ( r & c )- From ( 16 ) and ( 18 )> we S et 

(IJ\IJ Q ) F(m)=m\ J2J2^f Xf (R) (25) 

Reg, f ^ d f 

= ^Ex E x f (R)xU(R) = §^ (26) 



m!(n-l)! 
(m + n-1)!' 



(27) 



In getting from (25) to (26), the character %j m \{R) = 1 of the totally sym- 
metric representation (to) of the permutation group has been inserted, and 
the orthogonality relation of the characters has been used. To get to (27), 
d( m ) — X(m)(e) — 1 as well as rf( m ) = (to + n — l)!/(n — 1)!to! (see (15)) have 
been used. 

The result in (27) agrees with the result (10). It is one of the very few cases 
where the value of the integrals can be obtained group-theoretically without 
knowing the values of the individual £[c]'s. 



3.3 Primitive Integrals 

Integrals in which both symmetry groups Qi and Qj consist only of the identity 
e will be called primitive. This happens when all the indices i a in the set I 
assume distinct values, and all the indices % in the set J are also different. 
The corresponding diagrams have p dots each on both columns, and precisely 
one solid and one dotted lines connecting to each of the dots. The primitive 
diagrams for p < 3 are shown in Fig. 6, and the ones for p = 4, 5 are contained 
in Appendix C. 

Since Gi = Gj = e, it follows from (17) that N[c] — S(Q G c), where Q 
can be any element of S p . The primitive integrals (16) are simply £[c], one for 
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(a) p = 1 



Qe] §[(12)] 
(b) p = 2 



^[(12)(3)] §[(123)] 
(C)p = 3 

Figure 6: Primitive diagrams for (a) p = 1, (b) p = 2, and (c) p = 3. The 
identity element is everywhere denoted by e. 



each class c of S p . We may therefore use an element of each cycle structure to 
label the primitive integrals, as is done in Fig. 6. Diagrammatically, the cycle 
structure is translated into the loop structure of its diagram, as can be seen in 
Fig. 6. Using (18) along with (15) and the character tables found in Appendix B 
(note that df = X/( e ))> the primitive integrals forp < 3 can be easily computed, 
and the results are displayed in Table 1. The results for p = 4, 5 can also be 
found in Table 4 of Appendix C. 





c 


p = l 


p = 2 


p = 3 






i 


n 2 -2 


e 


n 


7l 2 -l 


n(n a -l)(ri 2 -4) 


(12) 




-1 


-1 




n(n 2 -l) 


(n 2 -l) (71^-4) 


(123) 






2 






n(n 2 — l)(ri 2 — 4) 



Table 1: Algebraic expressions for the primitive diagrams of p = 1, 2, 3. 
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Figure 7: Arbitrary stack diagram, ^(pi,p 2 , ■ ■ ■ ,Pt), of degree p — ^2l =1 Pi- 



3.4 Stack Integrals 

The stack diagrams (see Fig. 7) are direct integrals made up of disconnected 
lines of arbitrary multiplicities. As such, Q = e, and J differs from I only 
by relabelling. Using item 3 of Sec. 2.2, we may assume J = I. Hence stack 
integrals are integrals of the form (II\II). 

Let pi,p 2 , ■ ■ ■ ,pt be the multiplicities of the disconnected lines in a stack 
diagram. Then Gi = Gj = G = S Pl ® S P2 ® • • • ® S Pt , and N[c] is nonzero only 
when the class c is a direct product of the classes c, of the groups S Pi . In that 
case, 

t 

where ni(a) is the number of elements of S Pi in the class Cj. In other words, 

- nP /l , . ( 29 ) 

where the class a consists of ctj cycles of length j. Denoting the stack integral 
by S(pi,p 2 , ■ --,Pt), we get 

E(pi,p 2 ,...,p t ) = ( n pi! "^ ci ) ) ^ ci<g,c2<g "" <8)c *)- ( 3 °) 

ci,c 2l ... \i=l / 

All stack diagrams can be obtained by making the assignment fi — ► from 
each representation. In this way, we expect a same number of stack digrams 
as of primitive diagrams, or classes. Using the £ expressions obtained in the 
preceding subsection, the stack integrals for p < 3 can be computed to yield the 
expressions in Table 2. 
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5 (pi, ...,p P ) 

= n "( 3 ) = n(n+l)'(n+2) 

S ( 2 ) = n(n+l) S ( 2 ' X ) = (n-l)n(n+2) 
S(l,l)= 7 j^T 5(1, 1) 1) = „(,^-l)(r^-4) 

Table 2: Algebraic expressions for the stack diagrams of p = 1, 2, and 3. 









a a+1 . 


. 2a 






'i : 


b a 


a 










• c 2a 








■ C 2a Ci 


c a 



3.5 Special Double-Fan Integrals 

The index sets for the fully opened double-fan integrals [A a ] a [A b ] Q (Fig. 3(b) 
with N = 2a) are 



(31) 



Hence both Qj and Qj Q consist only of the identity e. As for Qi, it is given 
by S a ® S a , where the permutation groups S a act respectively on the b and a 
indices in /. The element Q maps Jq to J, i.e. Q = (1, a + l)(2, a + 2) ■ ■ ■ (a, 2a). 

The fully opened integral can be computed using (16), with N[c] given by 
(17) or (19). Thus, 

N ® = E E HQTRec) 
R&Si Teg. jQ 

= S(QRec) = ^(Q'ec), (32) 
Reg, Q> 

where the last sum is over every permutation Q' that sends all b indices in 
(31) to the positions labelled from a+1 to 2a, and similarly all a indices to 
the positions labelled from 1 to a. As a consequence, the allowed cycles of Q' 
must be of even length, and they can be specified by a sequence of nonnegative 
integers (k) = (fci&2 • ■ - ha), ki being the number of cycles of length 2i. The 
number of Q' with the class structure (k) that is related to c is given by 

w|c| " mS^- <33) 

In order to see how this is arrived at, consider an example where k\ — 2, k 2 = 2, 
and all other ki values are zero. Then Q' is of the form (ba)(ba)(baba)(baba), 
where the b and a letters should take the distinct index labels in (1, ... , a) and 
(a + 1, ... , 2a) respectively. Another Q' with the same cycle structure can thus 
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a 




1 

2 
3 


-1 
n(n 2 — 1) 
2 

(ri 2 -l)ri 2 (ri+2)(n+3) 
-6 


(ri-l)n 2 (n+l) 2 (ri+2)(ri+3)(ri+4)(ri+5) 



Table 3: Values of the monomial integrals [Aj] Q [Ab] Q for a = 1,2, 3. 




Figure 8: Non-orderly integrals of p — 3. 



be obtained by permuting individually all the a and b labels. This accounts 
for the numerator in (33). However, such permutations do not necessarily give 
distinct Q' elements. The cyclic nature of a cycle tells us that each cycle of 
length 2i will appear i times; this accouts for the i ki factor in the denominator. 
Moreover, no new Q' is obtained if we permute cycles of the same length; that 
accounts for the other factor fcj! in the denominator. 

We may now return to (16) to calculate the integral [A a ] Q [Ab] Q in terms of 
the primitive integrals £[c]. The result for the first few a values are listed in 
Table 3. 



3.6 Non-Orderly Integrals 

All integrals with degree p < 3 arc orderly. The non-orderly integrals of p = 3 
are shown in Fig. 8, and those related to them by the fan relation (8). The 
calculation of N[c] and the integral for each of them is discussed below. The 
integrals will be labelled by their figure, e.g., integral I(8a). 
The index sets for Fig. 8(a) are 



label \ 




2 


3 






b 


a 






c 


c 


Jq ) 


\d 


c 


c 



(34) 



They give rise to the symmetry groups Qi = {e, (12)} and Qj = Qj Q = {e, (23)}. 
Moreover, the element Q can be taken to be the identity clement. In order to 



15 



obtain the coefficients N[c] of equation (16), we need to compute QTR for all 
T e Gj Q and R e Gi- That triple product is QQj q Qi = {e, (12), (23), (132)}. 
As a result, 



/(8a) - Z(l, 1, 1) - £[e] + 2£[(12)(3)] + £[(123)] - 

In the same way, the index sets of Fig. 8(b) are 

' label ^ 
I 
J 

Jq 



(n 2 -l)(n + 2)' 




(35) 



(36) 



and hence the symmetry groups are Qj Q — {e} and Gi — {e, (23)}, and the 
exchange element is Q — (12). We thus obtain QGj q Gi = {(12), (123)}, from 
which 

f ' 8t '^'< 12 »< 3 " + ^" 123 "° ^-iK„ + 2) <37 » 

follows. 

Finally, for Fig. 8(c), the index sets are 

/label\ /l 2 3* 
/ 
J 

V Jq 

The relevant symmetry groups are Gj Q — {e, (13)} and Gi = {e, (23)}. With 
Q = (12), the set QGj q Gi is {(12), (13), (123), (132)}, and formula (16) gives: 

-2 




(38) 



/(8c) = 2(^[(12)(3)] + C[(123)]) = 



(n 2 - l)n(n + 2)' 



(39) 



The calculation of QGj q Gi is not that cumbersome for p = 3, but it gets 
worse pretty quickly as p increases. For example, let us look at some examples 
of p = 4. 

Let us first calculate Z(2, 1, 1) of Fig. 1(b), whose index sets are 



(40) 



Then Q = e, Gj Q = {e, (12), (34), (12)(34)} and£ 7 = {e, (12), (13), (23), (123), 
(132)}. Thus QGj q Gi = {e, (12), (13), (23), (123), (132), (12), e, (132), (123), 
(23), (13), (34), (12)(34), (143), (243), (1243), (1432), (12)(34), (34), (1432), 
(1243), (243), (143)}, hence 

Z(2, 1, 1) = 2^[e] + 8C[(12)(3)(4)] + 8£[(123)(4)] + 2^[(12)(34)] + 4^(1234)] 

= (41) 

(n-l)n(n + 2)(n + 3) 1 ; 



/ label \ 




2 


3 


4 \ 




b 


b 


b 


a 


J = 


d 


d 


c 


c 


V Jq ) 


\d 


d 


c 


e) 
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Figure 9: Sigma integral with all lines being of multiplicity one. 



Our last example is the E-integral, shown in Fig. 9. Its index sets are 



/ label 
I 
J 




1 

b 



e 



2 
b 
d 
d 



a a 



3 4\ 



d c 
d cj 



(42) 



V j q 



e 



Hence Q = e, Qj Q = {e, (23)} and Qi = {e, (12), (34), (12)(34)}. Multiplying 
such elements accordingly, the set {e, (12), (34), (12)(34), (23), (132), (234), 
(1342)} is obtained for QQj q Qi- Hence 

S = C[e]+3e[(12)(3)(4)] + 2e[(123)(4)]+e[(12)(34)]+e[(1234)] 

= H±l (43) 

(n- l)n 2 (n + 2)(n + 3)' 1 ' 

4 Comparison of the IM and the GTM 



We have discussed the computation of U (n) integrals (2) in two ways: the IM in 
Sec. 2, and the GTM in Sec. 3. Each of these two methods has its own merits, 
and drawbacks, and in a way they complement each other. The purpose of this 
section is to compare their relative stong and weak points. 

The IM is based solely on the unitarity condition (4), and the invariance 
of the Haar Measure (5). The method is simple because there is no need to 
know group theory. The conditions relate integrals of the same degree, and 
also integrals of degree p to integrals of degree p—1. Through these relations, 
specific integrals such as the fan integrals (10) and the Z-integrals (9) can be 
obtained, and general relation such as the single-fan relation (8) and the double- 
fan relation (12) can be worked out. 

The GTM has the advantage of being general, in the sense that all integrals 
can be computed using the formula (14) or (16). The price to pay is that we 
have to know the characters of the irreducible representations of the appropriate 
symmetric group, and a triple sum has to be carried out, which can prove to be 
very tedious for integrals of high degrees. Furthermore, unlike the IM, relations 
between integrals are hard to come by, so one must calculate the integrals one 
by one. There are however certain class of integrals, the orderly integrals, for 
which relations can be developed, and the quantity N[c] in (16) can be relatively 
easily computed. Then we merely have to know the primitive integrals £[c] in 
(16) to get the value of the orderly integral on hand. The stack integrals (30) 
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Figure 10: Unitarity sum relation involving the E diagram (second from the left). 
The sum is performed on the index indicated by an arrow. Using the fan relation 
(8), the unitary sum can be written as: + l)Z(2, 0, 2) + 2S = \Z{2, 0, 1). 



and the special opened double-fan integrals (33) are examples of this kind. The 
primitive integrals £[c] themselves must be calculated using (18). 

To summarize, the IM gives a huge number of relations but it is not easy 
to obtain the value of any specific integral. The GTM allows us to calculate 
any specific integral, albeit rather tedious at times, but it is difficult to ob- 
tain relations between integrals. In the next section, we shall discuss a hybrid 
method which makes use of the advantages of both methods. We shall use the 
general GTM formula to calculate a specific set of integrals, and then use the 
IM relations to obtain all the other integrals. 

In the rest of this section, we shall enlarge these general remarks about the 
IM and GTM, by using specific examples presented in the last two sections as 
concrete illustrations. 

The single-fan relation (8) can be obtained by both the IM and the GTM. 
However the double-fan relation (12) in its general form can be obtained only 
by the IM, because most of the integrals involved are not orderly, making it 
hard to derive relations using the GTM. Nevertheless, in special cases involving 
only orderly integrals, (22), the GTM can also be used to derive the relation. 

The Z-formula (9) is obtained using the IM, by a series of relations connect- 
ing it down to J dU = 1 . Since the ^-integrals are not orderly, it is hard to 
compute them using the GTM except at low degrees. The calculation of those 
by the GTM is shown in equations (35) and (41). 

However, since the values of the integrals in the IM are obtained only through 
relations, it may be relatively complicated to calculate just one specific integral. 
This is where the GTM is superior, because of the general formula (14) valid for 
any one integral. For example, it is easy to obtain the £ integral (43), assuming 
of course the £[c]'s to be already known. We can also obtain it using the IM, 
as we shall show below, but that involves a few steps because we must get it 
from relations. To see how that is done, look at Fig. 10, which is the unitarity 
relation applied to the dot of the first diagram indicated by an arrow. The first 
and third diagrams can be related to 2(2,0,2) by using the fan relation (8), 
and similarly the diagram on the right can be related to 2(2, 0, 1). Using the 
Z-formula (9), we then obtain 



S = \ Z(2,0,l) 



n — 3 



+ 1 2(2,0,2) 



2 
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n + l 



(44) 



(n- l)n 2 (n + 2)(n + 3)' 
the same as the result (43) obtained by the GTM. 



5 Hybrid Method 



Having understood the relative merits of the GTM and the IM, it is possible to 
combine their strengths into a more efficient hybrid calculational scheme. 

The strategy is to start with one or more integrals that can be computed 
by the GTM with relative ease. Generally speaking, such integrals are ordered. 
Once they are obtained, the many relations of the IM can be used to calculate 
other integrals from them. 

To illustrate this strategy, we will consider how the hybrid method can be 
used to calculate all double-fan integrals. 

By a double-fan integral, we mean any integral with two L-dots and any 
number of R-dots. Fig. 3(a) shows a closed (double-fan) integral (with the 
understanding that there are no extra dots or lines than those shown), and 
Fig. 3(b) shows a fully opened (double-fan) integral. We may also have partially 
opened (double-fan) integrals, in which every branch, namely, every R-dot with 
its connecting lines, can be regarded as a closed integral. See Fig. 11 for an 
example of a partially opened integral. 

As in Sec. 2.2, a fully opened integral is denoted by [A a ] a °- [A b ] ah [B a ]P« [B b \ f3b , 
and its corresponding closed integral is denoted by [a a A a +abAb+(3 a B a +f3bBb\. 
For a partially opened integral, we will denote it as a product of the closed 
integrals of each branch. See Fig. 11 for examples. 

Using (12) and (13), all double-fan integrals can be expressed as sums of 
fully opened integrals. Integrals of the form [A a ] Q [A;,] Q are given by (33) and 
Table 3, but we still have to know how to calculate a fully opened integral when 
[3i ^ 0. As shown in Appendix D, the IM allows us to relate them to those with 
Pi = 0, by using the following formula 



We close this section by showing how to use (12) and (13) to calculate the 
integrals in Fig. 11. 



[AX[Ab} a [B a f«[Bbf b 




(45) 



5.1 Fig. 11(a) 



There are two equivalent forms for this diagram. One is 



[A* + 2A b ][A a ] = {2[A a ][A b } 2 ) [A a ] = 2[A a ] 2 [A b } 2 



(46) 
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(a) (b) 

Figure 11: Partially opened double-fan integrals, (a) There are two equivalent 
forms for this graph: [A a + 2A b ] [A a ]and [A b + B a + B b ] [A a ] ■ (b) There are four 
equivalent forms for this graph: [A a + A b + B a ] [A a + A b ] , [2B a + B b ] [A a + A b ] , 
[Aa +A b + B a ] [B a + B b ] , [2B a + B b ] [B a + B b ] . 



and the other is 

[At, + B a + B b ] [A a ] = 4 ([A b ] [B a ] [B b ] ) [A*] - A[A a ] [A h ] [B a ] [B b ] , (47) 

where equations (12) and (13) have been used. The integral /(11a) is obtained 
by adding up (46) and (47). 

Using (45), we can express all fully opened integrals in the form [A a ] Q [Ab] Q . 
Applying to the present case, we get 

[Aa][Ab][B a ][Bb] - J^flMlM j^[A a nA b ]\ (48) 

Using Table 3, we finally obtain 

/(11a) = 2[A a ] 2 [A b ] 2 +A[Aa][Ab][B a ][Bb] 

-4 

(n 2 - l)n(n + 2)(n + 3)' 

5.2 Fig. 11(b) 

As shown in Fig. 11(b), /(lib) has four equivalent forms. For one branch, the 
factors are 



[Aa+Ab+Ba] = 4[A a }[A b }[B a }, 

[2B a + B b ] = 2[B a ] 2 [B b }; (49) 
and for the other branch, they are 

[A a +A b ] = [A a ][A b ], 

[Ba+Bb] = [B a ][B b ]. (50) 
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Hence 

7(116) - 4 [A a ] 2 [A b ] 2 [B a ] + 6 [A a ] [A b ][B a } 2 [B b ] + 2 [B a f[B b f. (51) 

We will now express each of the three monomial integrals in (51) in terms 
of [A Q ] Q [A 6 ] Q . First, with respect to (45), [A Q ] 2 [A 6 ] 2 [5 a ] is characterized by 
a = 2, (3 a = 1, and [3 b = 0. The vanishing of (3 b causes (45) to consist of the 
single term: 

[A^[A b Y[B^ = -^—[A a Y[A b ]\ (52) 
(n + 4) 

Second, [A a ] [A b ] [B a ] 2 [B b ] has a = 1, /3 a = 2, /3& = 1, and the sum in (45) gives: 
[A a U b ][B a f[B b] = {n + 2) l {n + 3) [A a ][M - (n + 2) 2 (n + 4) K] 2 [A] 2 - (53) 

Finally, [-B a ] 3 [S{,] 2 , having a = 0, /3 a = 3, (3 b = 2, can be expressed as 

1 6 
[Ba]3[Bb]2 = n2(n + l) 2 (n + 2) ~ n(n + 2) 2 (n + 3) [Aa] [A] + 

" \A a ] 2 [A b ? (54) 



(n+l)(n + 2)(n + 4) 



from equation (45). Using the fan relation, notice that [B a } 3 [B b } 2 can also be 
reduced to ^ ^Z(3, 0, 2). 

The expressions of [A a ][Aft] and [A a ] 2 [Aft] 2 in terms of n have already been 
determined in Example 1. The final answer is obtained by inserting (52)-(54) 
into (51). The result is: 

2 6(n — 2) 

[A a + A b + B a ] [A a + A b ] = n2{n + mn + 2) + n{n ; mr ! +3) {A a ][A b] + 

4(n2+2) WW 



(n + l)(n + 2)(n + 4) 
2(n 2 + 2n + 4) 



(n 2 - l)n 2 (n + 2)(n + 3)(n + 4) ' 
which can be verified using the plain group theoretical formula (14). 



6 Conclusion 

In this article, we have pursued the goal of finding an efficient method to calcu- 
late the monomial integral (1) or (2). We find that the IM discussed in Sec. 2 
is superior for deriving relations between integrals, but the GTM is able to give 
a formula to calculate any integral. The GTM formula involves a triple sum 
whose computation is often tedious and prone to mistakes. The sums simplify 
for orderly integrals, in which the invariant groups Qi and Qj Q are either dis- 
joint, or one is contained in the other. For non-orderly integrals, the hybrid 
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(c) (d) 

Figure 12: Diagrams that result from rotating away the pairs of lines (a) [A a ], 
(b) [At], (c) [B a ], and (d) [£?{,], from the compound pattern of Fig. 3. 

method is probably the most efficient. It uses the IM to relate them to some 
orderly integrals that can be calculated by the GTM with relative ease. 

This research is supported by the Natural Sciences and Engineering Research 
Council of Canada and by the Fonds de recherche sur la nature et les technologies 
of Quebec. 



Appendices 
A Derivation of Equations (12) and (13) 

To prove (12) and (13), we use the rotation technique discussed in item 5 of 
Section 2.1, and equation (7), to spin off from the R-dot of Fig. 3(a) a new 
R-dot attached to a pair of solid-dotted lines. Depending on whether the basic 
pattern of this new R-dot is [A a ], [Ab], [B a ], or [Bb], we get the graphs shown 
in Fig. 12(a), 12(b), 12(c), and 12(d) respectively. 

Repeating this spin-off operation over and over again on the R-dot still 
containing a compound pattern, eventually we come to a graph where every 
R-dot is given by a basic pattern. The resulting integrals are given in (12), 
corresponding to the decompositions of the compound pattern into a sum of all 
possible fully opened integrals obtained by spinning off. 
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The coefficient v(a' a a' b f3' a [3' b ) of these basic integrals is derived from a com- 
bination of three factors: 

1 . Each time that we spin off a basic pattern from a compound pattern with 
d pairs of solid-dotted lines, there is a factor 1/d arising from eq. (7) , by 
taking e = 1. Since we start from a compound pattern with N pair of 
lines, by the time we come to a fully-opened pattern we have accumulated 
a factor 1/AH. 

2. The N R-dots in the final pattern that is fully opened can be spinned off 
in a different sequential order. According to (7), they must be summed 
over. This gives rise to a factor N\/a' a \a' b \f3' a \/3' b \. 

3. At any time when we spin off a basic pattern, we can choose its single pair 
of solid-dotted lines in all possible ways. Eq. (7) says that we must sum 
over all these possibilities. The mulitipliciy factor is given by the number 
of permutations of these lines that lead back to the same basic integral. 
It is a factor of m a \ mt! n a l ml. 

Assembling these three factors, we get 

v(a a a b f3 a f3 b ) = , ngngn , (55) 
which is identical to (13). 

B Character Tables 

The character tables for p = 2, 3 are given here in the form used by M. Hamer- 
mcsh in [9]. The rows are labelled by the partitions that define the represen- 
tations, and the columns are labelled by the cycle structures that define the 
classes. The number of elements in each class, n(c), is written above the classes. 
The table for p = 1 is trivial, and it consists of the sole value 1. 

B.l p = 2 



Part.\Class 


1 1 

(I 2 ) (2) 


(2) 
(I 2 ) 


1 1 
1 -1 



B.2 p = 3 



Part.\Class 


1 3 2 
(I 3 ) (1,2) (3) 


(3) 
(2,1) 
(I 3 ) 


1 1 1 

2 0-1 
1 -1 1 



23 



C Primitive Diagrams for p = 4 and p = 5 



Using (14) and the character tables for 54 and S$ [9], the algebraic expressions 
for the primitive diagrams of Fig. 13 and Fig. 14 can be obtained, and they are 
given in Table 4. 



• *: P 




£[e] £[(12)(3)(4)] £[(123)(4)1 ^[(12)(34)1 £[(1234)] 



Figure 13: The p — 4 primitive diagrams. 



£M 



$[(12X3X4X5)] £[(123X4X5)] £[(12)(34)(5)] 



£[(1234)(5)] £[(123)(45)] £[(12345)] 

Figure 14: The p — 5 primitive diagrams. 
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p — 4 


p — 5 




n 4 -8n 2 +6 


n 4 -20n 2 +78 


e 


n 2 (n 2 -l)(ri 2 -4)(n 2 -9) 


n(n 2 -l)(n 2 -4)(ri 2 -9)(n 2 -16) 


(12)(3)(4) 


-1 


-(n 2 -2)(n 2 -12) 


n(n 2 -l)(ri 2 -9) 


n 2 (n 2 -l)(n 2 -4)(n 2 -9)(ri 2 -16) 


(123)(4) 


2n 2 -3 


2 


n 2 (n 2 -l)(ri 2 -4)(n 2 -9) 


n(n 2 -l)(ri 2 -4)(n 2 -16) 


(12)(34) 


n 2 +6 


n 2 -2 


n 2 (n 2 -l)(n 2 -4)(n 2 -9) 


n(n 2 -l)(ri 2 -4)(ri 2 -9)(ri 2 -16) 


(1234) 


-5 


-5n 2 +24 


n(n 2 -l)(n 2 -4)(n 2 -9) 


n 2 (n 2 -l)(ri 2 -4)(ri 2 -9)(ri 2 -16) 


(123)(45) 




-2(ri 2 + 12) 




n 2 (n 2 -l)(ri 2 -4)(n 2 -9)(n 2 -16) 


(12345) 
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n(n 2 -l)(ri 2 -4)(ri 2 -9)(ri 2 -16) 



Table 4: Algebraic expressions for the primitive diagrams of fourth and fifth 
degrees. In the p = 5 case, the elements from row two to row five should be 
written with the additional (5) one-cycle. 



D Derivation of Equation (45) 

We would like to express the general fully-opened integral [A a ] a [Ab] a [B a ]^ a [Bbf b 
in terms of the special ones of the form [A a ] a [Ab] a . The idea is to apply a uni- 
tarity sum on the [B a ] or [Bb] basic patterns to get rid of them. To get the 
final result we also need to apply the fan relation (8) or the double-fan relation 
(12) and (13). Our approach is to first determine how can [A a ] Q [A6] Q [_B a ]^ a 
be reduced to fully opened integrals involving only the [A a ] and [Ab] patterns. 
With such an information at hand, we will then try to reduce the more general 
[A a ] a [Ab] a [B a y 3 "[B b ] f3b integrals into the [A a ] a [A b ] a [B a f- integrals. 

Let us apply a unitarity sum on one of the [B a ] patterns in [A a ] a [Ab] a [B a ]P a : 

(n-(2a + f3 a - 1)) [A a ] a [A b ] a [B a f« + (& - 1) [A a ] a [A b ] a [B a f«- 2 [2B a ] 
+ a ( [A a ] a [A b ] a - 1 [A b + BaW'- 1 + [AX'^Aa + B a ][A b ] a [B a ] l3a ~ 1 ) 
= [A a ] a [Ab] a [B a f'- 1 (56) 

Using (8), [2B a ], in the second term above, can be rewritten as 2[B a ] 2 . Fur- 
thermore, equation (13) tells us that [Ab + B a ] and [A a + B a ] can respectively 
be rewritten as 2[A&][.B a ] and 2[A a ][_B a ], and the term in parentheses above 
simplifies to 4[A a ] a [Ab] a [B a ]P a . As a result, relation (56) reduces to: 

[A a ] a [Ab] a [B a Y 3a = n + 2a \ f3a _ l [A a ] a [Ab] a [B a f«-\ (57) 

Using the relation (57) recursively on its right-side, until no [B a ] remains, we 
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obtain: 



[A a ] a [A b } a [B a f° 

(n + 2a + f3 a - 1) (n + 2a + J a - 2) " ' (n + 2a) ^ ^ 
(n + 2a-l)! 



(n + 2a + f3 a - 1)1 



[A a ] a [A b } a , (58) 



and the first step of the work is completed. 

Assuming that p b < Pa, let us perform a unitarity sum on a [B b ] pattern in 
[A a nA b ] a [B a ]HBb] Pb : 

(n - (2a + a + p b - 1)) [A a ] a [Ab} a [B a }^[Bb} f:)b 
+ (p b - 1) [A a ] a [Ab] a [B a f«[Bbf b - 2 [2B b ] 

+ P a [AXlAbUBaf^lBa + B^B,}^ 1 

+ a ( [AXW-'iAt, + B b ][B a )^[B b ]^ 

+ [A a ] a -\A a + B b ] [A b ] a [B a f° [B^- 1 ) 
= iA a } a [A b } a [B a f°[B b f>-\ (59) 

Relation (8), or formula (13), again permit to make some simplifications, i.e. 
[2B h ] = 2[B h ] 2 , [B a +Bb] = [A a ][Ab} + [B a }[B b ], [A b +B b ] = 2[A b ][B h ], [A a +B b ] = 
2[A a ][B b ]. By making the proper substitutions in (59), the recursion equation, 

[A a ] a [Ab] a [B a f-[Bbf b = 

(n + 2aU-l) {[ ^ a[iUlaWfcW<fc " 1 

- P a {Aa\ a + 1 {Ab\ a + 1 [Baf a ~ 1 [Bbf i ~ 1 } , 

results. The above can be solved to give 

[Aa] a [A b ] a [Baf«{Bbf b = 

'Pa\ (Pb 



e=0 V 



e I \ e 



(n + 2a - 1 + 2e) • 
[A a } a+e [Ab} a+e [BaY 



{n + 2a-2 + e)\ r A + + Pa _ e 



(n + 2a + p b -l + e)\ 
which upon substitution of (58) yields the desired equation (45). 
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